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We consider a modified gravity theory, f{R) = R — a/R^ + hR"^, in the metric 
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Q> ■ formulation, which has been suggested to produce late time acceleration in the Uni- 
(N 

verse, whilst satisfying local fifth-force constraints. We investigate the parameter 

(N 

I range for this theory, considering the regimes of early and late-time acceleration, Big 

■ Bang Nucleosynthesis and fifth- force constraints. We conclude that it is difficult to 

. find a unique range of parameters for consistency of this theory. 

o 



Ph! I. INTRODUCTION 

Various cosmological observations suggest that the universe is pervaded by a new form 
^ \ of energy, dubbed dark energy, giving rise to accelerated expansion at the present epoch 
[1, 2, 3]. These observations are a challenge for fundamental physics, since a well-motivated 
candidate for dark energy has to be found. The usual candidates include scalar fields or 
extra dimensions. 

Recently, it has been suggested that instead of a new matter form, the cosmic accelerated 
expansion could be attributed to a modification of General Relativity itself, see e.g. [4, 5, 6, 
7, 8, 9, 10, 11, 12]. For reviews see e.g. [13, 14, 15]. The simplest models of modified gravity 
are theories in which the Ricci scalar R in the Einstein-Hilbert Lagrangian is replaced by 
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some (arbitrary) function of R. Examples of these models include inverse powers f{R) = 
or f{R) = R — aR~^. It has been shown that these models are equivalent to scalar tensor 
theories with vanishing Brans-Dicke parameter, thereby ruling out these models [16]. It 
has recently been argued that by adding terms like hR™ (m positive), the field can be 
made massive and hence short-ranged, allowing the theory to be made consistent with local 
constraints [10]. However, past work considered only the vacuum of the theory, without 
taking into account the coupling to matter. 

In this paper, we will consider theories with positive and negative curvature terms, taking 
the matter coupling into account. This is a non-trivial task and we have to make simplifying 
assumptions in order to proceed. The paper is organised as follows. In the next Section 
we formulate the theory both in the Jordan frame and in the Einstein frame. We will 
write down useful equations needed for later Sections. In Section III, local constraints for 
the choice f{R) = R + hR^ are considered, analytically and numerically for m = 2 and 
numerically for more general cases. In Section IV we discuss consequences for the model 
f{R) = R — aR~^ + bR"^. We will summarise our findings in Section V. 



II. JORDAN AND EINSTEIN FRAME FORMULATION 



The theory we are going to consider is specified by the action 

5jF = j d\^g l^-^fiR) + , (1) 

where R is the Ricci scalar and £m is the matter Lagrangian. From this action, the field 
equations can be easily derived: 

R,J'iR) - lg,J{R) + g,.Df'{R) - V,VJ'{R) = yT^„ (2) 

where /' = df /dR. The trace of this equation gives 

Spherically symmetric solutions of f{R) theories have been recently studied in [17, 18, 19]. 
In particular it was found that the Schwarzschild-de Sitter metric is an exact solution to 
the field equations for a large class of models. 
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Furthermore, it is well known that the theory can be rewritten as a scalar tensor the- 
ory (e.g. [10, 20]). To see this, introduce two auxiliary fields, A and B, and write the 
gravitational part of the full action as 

'^JF = ^ / ^'^v^ (BiR -A) + f{A)) . (4) 

Making the variation with respect to B gives A = R, whereas the variation with respect to 
A gives 

B = f'{A\ (5) 
allowing us to eliminate B and write the gravitational part of the action as 

^''' = ^1 ^'^v^ (/'(^)(^ - ^) + /(^)) • (6) 
Variation with respect to A gives 

f"{A){R-A) = 0. (7) 

This equation is solved by i? = A if /" is non-vanishing. 

To obtain a Einstein frame formulation of the theory, perform a conformal transformation 
dfiu = G'^giiu of the action (6) and choose 

a = -\nf{A). (8) 

The gravitational sector in the Einstein frame then reads (Einstein frame quantities are 
denoted with a tilde) 

Sef = ^ [ d^x^g (r-1 fVa)' - V{a)\ , (9) 



J ' ' ' " V ~ 2 
with 

where the relation between the field A and a is given by Eqn (8). 

A few comments are in order. Firstly, for the conformal transformation to exist, /' has 
to be non-zero and positive, as it can be seen from Eqn (8). We will always consider an 
f{R) such that this is the case. Secondly, in order to identify the field A with the Jordan 
frame Ricci scalar R, f" must be non-zero, otherwise the solution to Eqn (7) is not unique. 
However, we will assume that even at the (isolated) points in which /" vanishes, A = R is 
true. We will not distinguish between A and R for the remainder of this paper. 
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From these considerations it has been shown that the theory is equivalent to a scalar- 
tensor theory with vanishing Brans-Dicke parameter cjbd- This is in contrast with current 
observations, which demand cubd > 40000, if the field is long-ranged in the solar system. 
However, for 

f{R) = R-^ + bR"^ (11) 

with a and b positive and non-zero, it has been argued that the scalar degree of freedom a 
can be made massive, leading to a short-ranged force and thereby avoiding the confiict with 
local (i.e. solar-system) experiments [10]. 

Considering a fiat Robertson- Walker universe filled with dust, the Einstein-frame field 
equations give 

H' = ^ (P.n + P.) (12) 

and 

IdV (3k^ ,^ , 

In the last equation, (5 specifies the coupling between matter and the field a and is given 
by (5 = 1/2. The relation between j3 and the Brans-Dicke parameter cjbd is given by 
= 3/(20730 + 12), so that P = 1/2 corresponds to cubd = 0. The energy conservation 
equation for matter is given by 

+ 3Hp^ = 2(3&p^. (14) 

Since the theory was introduced to obtain an accelerating universe at the present epoch 
without to resort to an additional energy form, we fix 

a= (10-^2QgV)2{-+i), (15) 

following [10] and [11]. Note that, with a 0, Minkowski space is not a solution of the field 
equations. 

A. Properties of the Potential 

In this paper we are interested in the theory given by Eqn (11) for the function f{R), 

for which the potentials in the Jordan and Einstein frame have interesting properties. The 

effective potential for R (see Eqn (3)), fulfils 

dV}P _ 2f{R)-Rf'{R)-(3n^p 
dR 3f"{R) ^ ^ 
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where we have assumed a pressureless fluid with T = —p. 
The extrema of the potential are given by R, such that 

2f{R) - Rf'iR) - (3k^p = 0. (17) 

For gravitational stability (i.e. no instability for i?), we require that the extremum is a 
minimum, 

> (18) 



dR^ 

and using Eqn (10), we flnd 

< r < ^. (19) 

In regions of high curvature, such as during Big Bang Nucleosynthesis (BBN) or locally 
on Earth, we expect to be able to ignore the a term. In this regime, the minimum is specified 
by 

R^,^ + (2 - m)hRZr, - Pk^p = 0. (20) 

Since p is positive, a solution with positive R exists only when m < 2 or the magnitude of 
the second term is smaller than -Rmin- In the latter case we find that a good approximation 
for the minimum is given by -Rmin ~ (3hi^p. This, however, is not true for m < 2, but holds 
in general for R > hR^ . 

The effective mass at the minimum [R = Rmin) is given by 

2 ^ '^cS f (-Rmin) / (-Rmin) -Rmin 



/i 



dR^ 3/"(i?min) 

(3 - m)R^~'^ 2{m - 2)R f3K^p{2 - m)R'^-^ 

3m(m — 1)6 3(m — l)m 3(m — l)mb 

p2— m 
min 



3(m — l)mb 



(21) 



In the last two lines we have ignored the aR "-term. Note that the mass is dependent on 
the ambient matter density. 



III. LOCAL CONSTRAINTS FOR f{R) = R + bR^ 

In regions of high curvature, for example locally on Earth, one might expect the inverse 
curvature terms to be sub-dominant. In this section, we will study this regime of the theory 
and investigate the effect of an additional curvature term resulting in a local fifth force and 



6 



set a = 0. We note, however, that if we require only small corrections to Einstein gravity, 
we may assume that 

h < (22) 

to an order of magnitude estimate. This limit is consistent with the existence of a minimum, 
Eqn (20). For a range of 6, we discuss the strength and range of the resultant force, in order 
to test if this is a viable alternative theory of gravity. In doing so we will solve Eqn (3) and 
not the full field equation (2). Although limiting, this method is easier and provides insights 
into the predictions of the theory, taking matter couplings into account. It also allows us to 
study extended objects (like the Earth), instead of point particles. 

Fifth-force experiments have placed strong constraints on the strength of any deviation 
of the gravitational potential from that predicted by Einstein's theory of General Relativity. 
It is usual to assume a Yukawa potential for the fifth force, 

where a is a measure of the strength of the force relative to gravity, and A describes the 
range over which the force acts. The allowed strength of such a force is constrained for a 
large range of scales from 10^^ — lO^^m [21]. If this theory of modified gravity is to prove 
viable, it must satisfy the experimental constraints on all scales. 

For simplicity, we consider the simple system of the Earth sitting in the local solar system 
medium. We assume that the Earth has a radius of te = 6.7 x lO^mm and a constant density 
of 5.5g/cm^ (which corresponds to approximately k^Pe ~ 10~^^GeV^). The inter-planetary 
medium has an approximate density of 10"^^g/cm^ {i^'^Pss ~ 10"^^GeV^), which we also 
assume to be constant. We neglect the effect of the other solar-system bodies (in particular 
the sun), and the Earth's atmosphere (which we will justify later). 

In this simple system, it is possible to consider the force on satellites in orbit around the 
Earth. Geostationary satellites sit at approximately 6r^, corresponding to around lO^m. At 
this level, the strength of a fifth force is constrained to a < 10^*. 

If we consider a time- independent, radially symmetric solution, Eqn (3) becomes 

r/<«y_«|i^„^ (24) 



dr"^ r dr f" \dr J dR 
We shall only consider a subset of the theory and consider models in which m = |, 2, |, 3. 
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A. m = 2 



We begin by considering a specific case of our tlieory, where f{R) takes the form 

f\R) = R + bR^. (25) 

To check our results, we will consider solutions to both the Jordan frame R equation, and 
the equivalent Einstein frame equation for a. The range of h (as an order of magnitude) is 
given by 

h < lO^^GeV-^ 

as can be seen from Eqn (22). 



1. Jordan Frame R{r) 

In this simple case, the equation of motion for R{r), given by Eqn (24) becomes 

d^R 2 dR n ^ 

-rv + —r-f^^ = -f^ p 26 

dr^ r dr 

where p = l/VOb. We assume that R{r = oo) = -Rmin with R taking the value which 
minimises the effective potential, -Rmin = Pn'^p. Rq and Roo denote the minima in the test 
mass and background field respectively. 
A solution to Eqn (26) is given by 

„ , cosh ur „ sinh ur „ 

R = A + B—^ + i?min. (27) 



We consider two, solutions, exterior and interior to the test mass. In order to satisfy our 

dR 
dr 



equation, we require that ^ = at r = 0. The interior solution therefore takes the form 



fl,...= ''^-'^>'"""^'' +fio (28) 

pr 

where Ri is the value of the field at r = 0. 

Externally, we impose two boundary conditions: R{r = Vc) = Rc and R{r = oo) = -Roo, 
and the external solution becomes 

-Rext = 1- -Roo (29) 
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FIG. 1: Numerical solutions showing the evolution of R{r) for two choices of Vc- In the left-hand 
plot, is 10^ bigger than in the right-hand plot, so that ^ 1 and — > respectively. All 
other parameters are kept constant. The analytic solutions are shown with dotted lines and exactly 
follow the numerics. 



At the boundary between the test mass and the background, the two solutions and their 
derivatives must match. This condition yields the following expressions for Ri and Re'. 

Ri = (i?oo - Ro) (1 + /ir-e) e-'^'-^ + Ro (30) 
R. = i^?^±^ + ^-^^ [1 - (1 + /irj e--] (31) 

Given this form of a solution for R, we can see that there are two different regimes of 
behaviour. When fiVc » 1, Ri ~ Rc, and the field sits at its minimum inside the body. As 
the field approaches the boundary into the outside medium it begins to evolve very quickly, 
and soon settles into the exterior minimum. When fire <^ 1, -Ri ~ -Roo and the field remains 
at a value close to the external minimum. Two such cases have been modelled numerically 
and compared to the analytics above, see Figure (1). The reader may note the similarity of 
these solutions to the thick and thin shelled regimes described in the chameleon model [22]. 
We define a parameter, 6 to differentiate between thick (5 —* 0) and thin (5 — * 1) regimes: 

Sk = (32) 



2. Einstein Frame (T(r) 

Rather than formulating the theory in the Jordan frame, we can consider the equations 
for the scalar field a{r) in the Einstein frame. We will see that this has some significant 
advantages later on. 
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The equation governing the evolution of a in the Einstein frame is given by 

^^2^ l<>^^ (33) 

dr-^ r dr 6 da 

where 

dV^ ^ Rf'-2f + PK^p ^ -R + Pk^p 
da {f'f (l + 26i?)2 ^ 

This assumes that the conformal transformation relates i? to cr by 

e"'^ = f\R) = 1 + 2bR. (35) 

We have already stated that, in order to study only small deviations to Einstein gravity, 
we require bR < 1. It is therefore appropriate to assume /' is close to unity and we can 
expand e"'^ ~ 1 — a. Hence 



a ^ -2bR. (36) 



In this regime, Eqn (33) becomes 



+ _ ^ _ = (37) 

dr"^ r dr 6b 3 

We again require that ^ = at r = 0, and that a sits in the minimum of its effective 
potential far away from the test body. The solution to this system is equivalent to Eqn (27), 
taking a = —2bR, as it can also be seen from using Eqn (36) in Eqn (37) to get (26). Similar 
to Eqn (32), we can define: 

6. = (38) 

3. Fifth Force 

In order to calculate the modified gravitational force in the physical frame, we begin by 
considering the force generated by the equivalent scalar field in the Einstein frame. 

We can calculate the strength of force which would be mediated by a in the Einstein 
frame, 



Fe = -l3MVr(y (39) 
Ve = (3Ma. (40) 
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The force measured in the Jordan frame is related to the Einstein frame force via the 
conformal transformation 

Fj = c-^Fe = f'FE (41) 

Note once again that the assumption of perturbed Einstein gravity results in /' 1, and 
therefore 

Fj ^ Fe, (42) 

ignoring higher order corrections. 

The external solution for R is given by, 

{Rc — Rod) rcC'^^^"^'^"^ 

Rcxt = 1- Roo, (43) 

r 

and so using Eqn (36) 

a{r)^-2b{R,-R^)r,^^^^y^ + aoo. (44) 

Note this is the same equation we might have obtained directly from Eqn (37). 

Ignoring any additional constants, the potential of the field can be immediately seen to 

be 

^ ' ^ Ro rl 87iMj,ir ^ ' 

where gj^^^ = Using the definition of a Yukawa potential (Eqn (23)) we find 



32j,(-^c Roo) 6 '^^ 



and 

X = - = V6b. (47) 
/i 

We now consider the two limiting cases of the theory. Firstly, for fire <^ 1 {Sr 0) we 
have e'^"''" 1 and Eqn (30) can be approximated by 

R, = R^ + ifir,)\Ro - R^)/3. (48) 

Substituting this result into our expression for a, and assuming that Rq ^ Roo gives 

a = Al3%^i = I. (49) 
6 
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Note that while this result is independent of b, by assuming the conformal transformation, 
/" 7^ and 6 7^ is implicit. 

The second limiting case occurs when fiVc ^ 1 {Sr — ^ 1) and Eqn (30) becomes 

Rc = (50) 

and for Rq > i?oo, a is now given by 

a = 6f3%—. (51) 
In this case a depends on b and Vc- For a given experimental radius, r^, there is a minimum 

2 

value of a which will occur for b = The minimum value of a for a given radius is 

«min = T7T = 0.46. (52) 
16 

2 

Note, however, this minimum may not be seen, since b = ^ may not occur within the regime 
for b in which fire ^ 1. 

It can therefore be noted that the minimum force for the complete range of b will be 
given by a = |, when Sr — * 0. For the upper limit of 6 = lO^'^GeV"^, A ~ lO^^m or 25rE- 

For large b, such a force would be ruled out by experiment. As b is decreased, a increases, 
but since the range decreases, the force will eventually become undetectable. As the range is 
reduced, however, it becomes necessary to consider a smaller scale experiment. For example, 
laboratory experiments measuring distances less than 1cm consider a test mass in a vacuum. 
This alters the experimental setup considered here and is beyond the scope of this paper^. 



B. General m 



We continue our investigation by considering more generalised formulations of the theory, 
with 

f{R) = R + bR^ (53) 

where m is a positive constant. 

Analytical solutions for the general case cannot be found for Eqn (24). It is therefore 
necessary to solve specific cases numerically. The "width" of the potential minimum is 

^ In addition, in reducing b considerably, we would have to include the a/ R term in the theory, which makes 
the analysis of local constraints considerably harder. 
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described using the second derivative of the potential, as in Eqn (21). As the density 
decreases, the mass increases and the "width" of the minimum decreases. We may consider 
the scalar field to be moving along an inverted potential, as in bubble nucleation [23] and, 
in this notation, the minimum should be thought of as a maximum (see [22] for a good 
description). Therefore, as the density decreases, the maximum increases in height but 
decreases in width. Numerically, this provides a challenge: an extreme fine-tuning of initial 
conditions is required in order to sit the field on such a maximum (or conversely, to place 
the field in its minimum). Even with high-precision calculations, solutions are difficult to 
obtain. It is numerically easier to consider the evolution of the scalar field a directly and to 
observe the solutions as the external density is decreased. 
The equivalent equation for a is given by 

i^ 2J^_l^^, (54) 

dr^ r ar 6 da 

If we restrict ourselves to the regime where 

R > mhR^ (55) 

we can make the approximations 

/' ^ 1 a ~ -mhB!^-^ (56) 
and so Eqn (54) can be approximated by 

dr"^ r dr ~'~ 3 rah ) 3 ^ ^ 

Note that a is negative and hence the third term is real. In general it is not possible to find 
analytic solutions to this equation, but numerical solutions are possible. As for -R(r), we 
shall denote the minima of a internally and externally as cxo and cToo respectively. At r = 0, 
a = Gi. 



1. Fifth Force 

Heuristically, with knowledge from the m = 2 case, it is possible to estimate how the 
strength of the fifth force will change with b. When 6^ 0, we might expect the strength to 
take a value independent of b. As b is decreased, 6„ increases until 6„ = 1. As this happens. 
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the difference between minima (o"o — (Too) becomes larger, while the transition length becomes 
shorter. This implies that the resultant force becomes stronger, but the interaction range 
decreases. 

From this argument, the most detectable force will occur for the maximum allowed value 
of b, when 6cr —>■ 0. We will initially consider this case, in order to calculate the strength 
and range of the force. 

We will consider four distinct regimes of Eqn (57); one internal solution (r < r^), and 
three external solutions (r > r^), which we shall denote by I, II, III and IV respectively. These 
correspond to neglecting combinations of terms in Eqn (57). 

Let us consider the internal solution first. If (5o- ^ 0, the system is heavily damped and 
Eqn (57) can be approximated by 

d'^a 2 da l3iPp 



dr'^ r dr 3 

A solution can be found 



0. (58) 



Mr) = ^ + a. (59) 

where we have used the condition that ^ = when r = 0. Note once again that a is 

negative and hence the magnitude of a decreases. Outside r = r^, the density drops and the 

system can initially be re-expressed by 

d'^a 2 da 
dr"^ r dr 

which has a trivial solution 



0. (60) 



aij{r) = ^ + C. (61) 
Matching the derivatives at r = rc of aj and an leads to an expression for A: 

A = (62) 



Assuming that A is small and using Eqns (23) and (40), A can be related to a: 

QPA 2 , 1 
« = ? = oP = (63) 



Hence, in this limit, the strength of the fifth force is seen to be independent of m or b. This 
value for a is the same limit as seen in Eqn (49) for m = 2. 

The third and fourth regimes occur when the potential terms (last two terms on the LHS 
of Eqn (57)) become important. Firstly, consider that both terms are of similar magnitude 
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m b r/r^ 



3 ^q25 ;lq15 

3 10^°^ 10^ 



TABLE I: Minimum decay radius for the Earth with varying m given by Eqn (65). Values for b 
are taken close to the upper limit and correspond to the numeric parameters in Figure (2). 

(regime III). We can approximate the radius of decay from a = | by the interaction range 
of (7, A = 1/yU, where ^ is given by Eqn (21). As the density decreases, the mass increases 
for m < 2 and decreases for m > 2. The interaction range, A, does the converse. 

Finally, consider the case in which the density is low enough that the final term can be 
ignored: 



Then the solution starts to decay when then the second and third terms are of the same 
magnitude. We make the assumption that cr(r) follows an up until this time, where C is 
related to the density and is negligible. It is a simple rearrangement to show that the two 
terms roughly equate when 



This radius sets the length of decay. 

Note that regimes III and IV are exclusive: either the density term contributes or not. 
If the density term contributes, the solution follows regime III and the decay scale is set 
by the mass term. If the density is below a threshold, regime IV, it can no longer affect 
the decay length, which is then set by Eqn (65). Hence this latter radius sets the minimum 
decay length and we shall call it Tmin. Between and rmin, we expect the strength of the 
force to be i. Table I gives the minimum radius of decay for varying m in the Earth system 
when b is close to the upper limit. 

In order to verify this behaviour, we ran numerical simulations within the Einstein frame, 
using Eqn (57). We ran the simulation for varying m, once again using values for b close to 




(64) 




(65) 
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the maximum range. The code exphcitly solved for a{r) and da{r)/dr allowing for direct 
calculation of the field's force, F„. 

As in Section III A, we chose the Earth system. For each choice of m, the external density 
was varied to observe the effect on the decay length. A physical situation is one in which 



The results are shown graphically in Figure (2). Numerically, a binary search was un- 
dertaken in order to find the value which meant the external solution reached a^o- No 
assumption was made as to the value of cTj. However, in all cases, Ci <^ ctq, leading to models 
in the thickshelled regime (as expected for large b values). The analytics described above 
should be a good approximation to the numerical solution. 

Regimes I and II given analytically above are easily seen for all choices of m. From Table 
I, it is justified that we do not see a decay for m = 3/2 for the range of r in our simulation. 
Only regimes I and II are visible and the force remains constant outside r = r^. For m = 2, 
the decay radius occurs around r ^ 8000, which is approximately equal to A = from 
Section III A 3. Note that the decay radius is independent of the density, as expected from 
the mass term given in Section IIIAl. Regimes III and IV are explicitly seen for m = 5/2 
and m = 3. As the density decreases, the decay length also decreases, until the density 
finally has no more effect. The minimum decay length matches Eqn (65) within an order of 
magnitude. 

Previously we argued heuristically that the force is a maximum when b takes its maximum 
value. In order to test this, numerical solutions were found as b decreased away from this 
regime. Due to numerical difficulty, it was not possible to simulate the complete range of b 
for each choice of model m. Instead b was increased over an order of magnitude to test the 
effect on a. 

In the general m model, it is not possible to show analytically that the force has Yukawa 
form. Indeed, the evolution of a does not fit perfectly to a Yukawa field, but an approxi- 
mation can be made. We will describe the force given by an effective strength, aes, which 
varies with radius. We define 



where a{r) is calculated numerically. A pure Yukawa field will yield a constant aes- The 
numerics directly calculate the first and second field derivatives with respect to r, a' and a" 



the external density is that of the solar system, H?pss ~ 10 ''^GeV^. 




(66) 
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FIG. 2: The ratio of forces due to the scalar field, a and Newtonian potential, shown for different 



m. 



and the effective mass, fi{r), is given by 



1 la'' lV-4+(^0' 

/i r = -— ± 

r 2 a' 2 r 



(67) 



where we take the positive root. It is straightforward to rearrange Eqn (66) to determine 
acfr(r) numerically. The results are shown graphically in Figure (3) for m = |, 2, |, 3. 



When b takes a maximum value (allowed by the condition for a minimum), 6„ ^ and 



a(r) 



The tails at high r values correspond to the field evolving from the minimum 



due to extreme fine-tuning of initial conditions. At this point, the field no longer follows a 
Yukawa potential. As b is decreased from its maximum, the effective field strength increases, 
as expected. For m = 2, acs is given to high precision by Eqn (46). For cases approaching 
the thinshell regime, the field decays quickly, reaching a^o in shorter and shorter lengths. 
The numerical calculation of /i(r) agrees very accurately with the mass given in Eqn (21). 
Since the interaction range is the inverse of it is directly proportional to Vb and shortens 
with decreasing b. 

The ranges of A for the numerical results are given with Figure (3). The minimum 
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FIG. 3: Numerical calculation of the effective field strength, Ocff as a function of radius. A constant 
a(r) denotes a pure Yukawa regime. Re = IQ-^^GeV^, Roo = IQ-^^GeV^. For the numerics 
shown: X/rs = 24 - 7600(m = |), A/r^ = IQ-^ _ 700(m = 2), X/rs = W'"^ - 70(m = |), 
X/rE = 10-1 - 3800(m = 3). 



interaction range shown is quoted as 10~^r£;, which corresponds to 69km. As b decreases, 
it may become possible to hide the field, since even though the field strength increases, the 
interaction range decreases. 

Note that until this point we have neglected the Earth's atmosphere. If the field always 
stays close to Rss (thick regime), the presence of the atmosphere makes little effect on the 
field and we may neglect the atmosphere in this regime. We have checked this assumption 
numerically in the case m = 2. If i?j ^ Rq (thin regime), the field rolls quickly from 
its minimum within the Earth to that of the Solar System. The transition length can be 
the same magnitude of the atmospheric depth and we would expect the field to sit at its 
minimum, both within the Earth and the atmosphere. The fifth force constraints upon this 
system would therefore arise from both the Earth-atmosphere boundary and the atmosphere- 
space boundary. Given that we have shown that a > 1/6, independent of the values of Rq 
and Roo, we can assume that our results will still hold. 
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Once again, as in the m = 2 case, when A becomes very small, local experiments testing 
a fifth force require a simulation of test cases in vacua, which is outside the scope of this 
paper. 



As mentioned already earlier, it has been suggested that a modified gravity theory with 
both positive and negative curvature terms can provide a mechanism for early and late time 
infiation, while hiding a local fifth force arising from the additional scalar degree of freedom. 
In this section, we consider the effect of these terms, both locally and cosmologically, and 
investigate the allowed parameter range of b. 



At early times in the universe, with extreme high curvature, the theory tends to the limit 
R + bR"^. It has been postulated that the higher curvature terms may drive early time 
infiation (c.f. R^ infiation [24]). It is possible to estimate the required value for b, such 
that fiuctuations in the Cosmic Microwave Background assume the observed amplitude and 
spectral index. For m = 2, b ^ IQ-^^GeV^ [25]. Of course, it should be remembered 
that, whilst it is advantageous that this theory may lead to early and late time infiation, 
it is certainly not a requirement. Therefore b need not be constrained by such results. We 
merely mention it for completeness. 



In regions of high curvature, f{R) can be approximated hjR + bR"^ as in Section III. 
As the curvature decreases, the a term becomes more significant, which eventually results 
in the loss of exterior minima. We have confirmed numerically that, provided an exterior 
minimum exists, the evolution of R{r) follows exactly the results obtained in Section III, 
even when a 7^ 0. 

In Section II, it was stated that /" > for stability and locally 



IV. DARK ENERGY AND f{R) = R- a/K'' + bR 
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A. Early Time Inflation 



B. Locally 



a 




GeV < 1. 



(68) 



b R ("+™) 
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For the Earth, the field assumes a minimum Re ~ 10 ^^GeV^ and the existence of minimum 
^ leads to an upper limit on b (Eqn (22)). Hence 

lQ(53m-3in-84) QeV'^^""^) < b < (IQ-^^^) GeV^^^""^) (69) 

within orders of magnitude. In the solar system, assuming R sits in the minimum, Rss ~ 
lO^'^^GeV^ and therefore 

^Q(78m-6n-84) GeV'^^""^) < b < (10"^^) ^'^""^ GeV'^^"™) (70) 

It is straightforward to show that, for the cases of m and n we consider, there is no overlap 
between the upper limit of Eqn (69) and the lower limit of Eqn (70). Therefore, for a given 
b, we cannot simultaneously find minima within the Solar System and the Earth. 

In the absence of an external minimum, we are unable to specify the value that the field 
takes at infinity. In this case, we can either assume an external field value, or say nothing 
regarding the evolution of R{r). 

We shall consider two exclusive regimes; we take b such that (a) a minimum exists in the 
solar system and (b) a minimum exists in the Earth. 

(a) For a minimum to exist in the solar system, the value of b is high. For this range of b, 
the inverse curvature plays no significant role and the results of Section III hold. Due 
to the large value of b, the theory is restricted to a thickshell regime, which implies 
R ~ -Roo always. In this case, we expect a = ^ within a minimum radius given by 
Table (I). 

(b) We consider a value of b so that a minimum exists on Earth and locally the theory 
takes the limit R + bR^. For a given b value, there will be an exterior field cut-off, -Ra 
for the existence of minima. For R above this cut-off, the theory is well described by 
the theory discussed in Section III. Again, we can expect a > |. 

Below the field cut-off, no exterior minimum exists and we are unable to predict R^o- If 
we assume, however, that the field sits at a value -Rext (N.B. not given by the minimum) 
, we can make some predictions about the theory. Firstly, if i?ext > -Ra then we have 



^ For m < 2, the upper limit is not imposed due to existence of minima, but we will see later that Big Bang 
Nucleosynthesis constraints force this limit. 
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the same case as in the paragraph above. Secondly, if i?cxt < Ra, the theory tends to 
R — a/R^ at large distances from the Earth and this is ruled out experimentally [16]. 

For a b value outside of these two regimes, we are unable to say where the field will sit at 
any given radius and we cannot calculate the strength of the force. 



C. BBN Constraints on b 



One of the basic pillars of modern cosmology is Big Bang Nucleosynthesis (BBN) and 
the modified gravity theory we consider has clearly to reproduce the observed light element 
abundances. 

Stringent limits on the abundances of ^He and ^Li lead to a constraint on the deviation 
of the Hubble parameter during BBN from standard cosmology. Generally, a non-standard 



expansion rate can be parameterised by an expansion rate factor S 



H 

Hgr 



where H is the 



expansion factor in the physical (Jordan) frame and Hgr is the expansion rate in Einstein's 
theory for the same matter content. This in turn can be related to the number of extra 
relativistic neutrino species, AA''^ [26], by 



S 



1 



7AiV, 
43 



1/2 



(71) 



It then follows that, for a deviation from General Relativity 



1 ^ 7AiV,/4 
52 10.75 + 7AiV^/4' 



In the Jordan frame, the variation of the Hubble parameter can be related to the conformal 
factor A = exp(a/2) [27] 



^BBN 



(72) 



Using Eqn (8) we find 



0.86 < < 1.19 



(73) 



where we have used jAA'^^l < 1 [26, 28]. The method above has been used to look for 
deviations of BBN in scalar-tensor theories. 
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During the BBN epoch, we expect positive curvature terms to dominate, due to the high 
curvature'^. At the present time, however, the observation of dark energy leads us to assume 
that the term involving inverse powers of R is beginning to dominate. These assumptions 
lead to the following constraint: 

1 + bmRT~^ 



1 + ani?-^"+'^ 



1, (74) 



where we have assumed that during BBN the a/R term is negligible and that the bR^ 
term is negligible for the cosmological evolution at the present epoch. If we require that 
acceleration today arises due to the inverse curvature term, the terms in the denominator 
must be of order unity. Hence, the b term must be order unity for the BBN constraint to 
hold. This leads to a fine-tuning of the parameter b: 

b ~ Rl-^^ (75) 

Note that this is a product of the coincidence problem of dark energy; the fact that the 
a term just starts to dominate today implies that the R^ term stops dominating at the 
BBN epoch. Since, we expect R to decrease since the BBN epoch, we can assume that such 
positive curvature terms will be sub-dominant after this time {bR"^~^ < 1). 

Therefore the evolution of the universe must closely resembles that produced by Einstein 
gravity at times from BBN until today. We can therefore approximate R = jSn^p, which is 
the condition for the field to be sat close to its minimum (see Eqn (20)). With this argument, 
we can estimate -Rbbn ~ '^^Po(l + -^bbn)^- Hence, 

b ~ (10^^^ GeV')^-"*. (76) 

Note that this value lies just inside the range for the existence of a minimum on the Earth. 



D. Present Day Acceleration 

As argued in the previous section, we should assume for cosmological evolution that the 
field sits close to its minimum (and f{R) ~ R) until today when inverse curvature terms 
begin to play a role. It is interesting to consider the constraint on b coming from this 



Deviating too much from this might result in a non-standard cosmological evolution and is ruled out [29] . 
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requirement. Note that, at recent times, we expect the field to move out of its minimum in 
order to induce acceleration. 

We numerically solve Eqn (17) in order to find R, using the full f{R) = R — a/R^ + bR^ 
theory. We always take the positive, non-zero root, since R > 0. Using this value, we 
solve Eqn (19) to find the limits on b. The current cosmological energy density is po = 
1.3 X 10~^''GeV^ and we assume a to take the value in Eqn (15). The limits for the cases 
we have considered are: 

(77) 
(78) 
(79) 
(80) 

Note that for m < 2, the upper limit in Eqn (19) is always satisfied, while for m > 2, the 
limits tightly constrain b. We also find that these limits are largely independent of n, over 
the range n = ^, 2. 



m = 


3/2 


b > lO^^GeV-^ 


m = 


2 


b > lO^^GeV-2 


m = 


5/2 


b ~ lO^^^GeV"^ 


m = 


3 


b ~ lO^^^GeV"^ 



V. COMPARISON OF REGIMES 



In this paper, we have looked at the many various constraints that can be placed on 
a modified gravity theory, f{R) = R — a/ R^ + bR"^. These constraints include local fifth 
force experiments. Big Bang Nucleosynthesis and cosmological evolution. In addition, we 
compare to early time infiationary results. 

For densities corresponding to the Earth, solar system and cosmology {pE,PsSyPo)y there 
is no value of b that simultaneously satisfies the existence of a minimum in all media. Indeed, 
b values exist for which there is an absence of minima in all these regimes. For these values, 
it is not possible to comment on the evolution of R either locally or cosmologically. 

We may choose b such that a minimum exists for one or other of these densities. These 
cases were discussed individually in Section IV and summarised schematically in Figure (4). 
We now discuss these regimes with reference to each other to judge the consistency of the 
model. 

Firstly, we can consider large values of b such that a cosmological minimum exists until 
very recently. At high curvature, for example on Earth, such large b values result in a 
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bR"^ theory. This theory has been considered previously [30], where it was concluded that 
experimental observations require 0.25 < m < 1. 

Secondly, when a minimum exists in the solar system, the regime is one of a thickshell, 
which can be ruled out for the models considered due to a = | and the minimum decay 
radius in Table I. 

Thirdly, the complex consideration of minima on Earth was discussed in Section IV. The 
most viable values for b exist at the lower end of the Earth bounds, when A is small enough to 
hide the scalar degree of freedom. However, it should be noted that the arguments presented 
in this paper rely upon a specific (large scale) experimental setup. To constrain the theory 
on scales A < Im would require a small scale experimental setup to be considered, (i.e. a 
system of test-mass, vacuum flask, and exterior medium). Furthermore, we can compare 
these lower b values with the BBN result in Eqn (76). The two constraints are inconsistent. 

It is of course possible that a viable cosmology could be produced without assuming that 
the field sits in the minimum of its effective potential. However numerically solving such 
a system is an extremely difficult problem, and we have been unable to produce a realistic 
cosmology. For a specific model (m = 2,n = 1), it was noted that if the field starts with 
large R (at early times), the system was not found to evolve into its local minimum [29], see 
however [31]. Note that the BBN result in Section IV C specifies that the field must exist 
close to its minimum from that epoch until recently. This is inconsistent with a cosmology 
in which the field does not settle to its minimum. 

Further to this discussion, if the model is required to produce early time inflation, none 
of the regimes above are consistent with the result found in [25]. 

VI. CONCLUSIONS 

In this paper, we have investigated the parameter range of the modified theory of gravity, 
f{R) = R- a/R"" + 

For the specific case in which a = 0, we modelled fifth force experiments using the local 
Earth system and calculated the strength and range of such a force. We find that, assuming 
a Yukawa form for the potential, the force has a > | and = ~ where /i is approximated by 
Eqn (21). This restricts the allowed parameter space for b such that the interaction range 
is smaller than experimentally allowed ranges. 
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FIG. 4: Schematic diagram showing the relations between regimes for 6, A and a. 

For the general model (a 7^ 0), we assume that a/i?" leads to late time acceleration, 
setting a for this to be the case. We considered the allowed range of b required by the 
following regimes: early and late-time acceleration. Big Bang Nucleosynthesis, existence 
of minima and fifth-force constraints. Our results show that it is very difficult to find a 
consistent parameter range to satisfy two or more regimes. 

The tightest constraint seems to come from BBN. If the curvature values are known 
today and during the BBN epoch, the parameter b would be constrained to a high degree. 
Comparison of this result with fifth force constraints can ultimately decide the viability of 
this theory. Of course, more complicated cases for f{R) might be considered, which might 
evade the afore-mentioned problems, see e.g. the discussion in [32]. 

In future, one should also verify our results by solving the full field equations in the 
presence of an extended object, like the Earth, and study the stability of the solutions, 
similar to [33, 34]. 
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